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Abstract 

We investigate the effect of antikaon condensed matter on bulk viscosity in rotating neutron 
stars. We use relativistic field theoretical models to construct the equation of state of neutron 
stars with the condensate, where the phase transition from nucleonic to K~ condensed phase 
is assumed to be of first order. We calculate the coefficient of bulk viscosity due to the 
non-leptonic weak interaction n ^ p + K~ . The influence of antikaon bulk viscosity on the 
gravitational radiation reaction driven instability in the r-modes is investigated. We compare 
our results with the previously studied non-leptonic weak interaction n-\-p ^ p + K involving 
hyperons on the damping of the r-mode oscillations. 

We find that the bulk viscosity coefficient due to the non-leptonic weak process involving the 
condensate is suppressed by several orders of magnitude in comparison with the non-supcrfluid 
hyperon bulk viscosity coefficient. Consequently, the antikaon bulk viscosity may not be able 
to damp the r-mode instability, while hyperon bulk viscosity can effectively suppress r-mode 
oscillations at low temperatures. Hence neutron stars containing condensate in their core 
could be possible sources of gravitational waves. 

1 Introduction 

Immense information about the internal composition of neutron stars can be obtained through 
the study of unstable modes of oscillations associated with rotating neutron stars. Among the 
various possible instabilities, inertial r-modes restored by Coriolis force are particularly interest- 
ing as they are thought to play an important role in regulating the spins of newly-born neutron 
stars as well as old, accreting neutron stars in low mass X-ray binaries (LMXBs). If the r-mode 
is unstable, a rapidly rotating neutron star could emit a significant fraction of its rotational en- 
ergy and angular momentum as gravitational waves, which could be detectable by the upcoming 
generation of gravitational wave detectors. 

Bulk viscosity in neutron star is caused by the energy dissipation due to non-leptonic weak 
interaction processes in pulsating dense matter. It was argued that the r-mode instability could 
be effectively suppressed by bulk viscosity due to exotic matter in neutron star interior. Neutron 
star matter spans a wide range of densities, from the density of iron nucleus at the surface of the 
star to several times normal nuclear matter density in the core. Since the chemical potentials of 
nucleons and leptons increase rapidly with density in the neutron star core, different exotic forms 
of matter with large strangeness fraction such as hyperons, Bose-Einstein condensates of antikaons 
or deconfined quarks may appear there. The coefficient of bulk viscosity due to non-leptonic weak 
processes involving hyperons was calculated by several authors [H EJ [3l HJ [5] . Also, the impact 
of bulk viscosity due to unpaired and paired quark matter on the r-mode instability had been 
investigated extensively [SI EJ [7] . 

In this paper, we study the effect of the presence of antikaon condensates on the equation 
of state (EoS), bulk viscosity and the corresponding damping timescale. We compare the bulk 
viscosities generated by the non-leptonic weak interaction n ^ p + K~ involving antikaons with 
that due to the non-leptonic process n + p ^ p + A. involving hyperons as previously obtained. 
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and investigate their role in damping of r-mode oscillations in rotating neutron stars. 

This paper is organized in the following way. In Sec. II, we describe the model to calculate 
equation of state, bulk viscosity coefficient and the corresponding timescale. The parameters of 
the model are listed in Sec. III. The results of our calculations are discussed in Sec. IV, and the 
summary and conclusions are given in Sec. V. 



2 Theoretical Model 

We assume a first order phase transition from nuclear to antikaon condensed matter. A relativis- 
tic field theoretical model is adopted to describe the /3-equilibrated and charge neutral matter in 
both the phases. 



2.1 Hadronic Phase 

The constituents of hadronic phase are octet baryons, electrons and muons. In this model, baryon- 
baryon interaction is mediated by the exchange of scalar and vector mesons. For hyperon-hyperon 
interaction, two additional strange mesons, scalar /o (denoted by a*) and vector <j) are incorpo- 
rated. The Lagrangian density for the hadronic phase is given by 

J~-B = (hfid^^ -1713 + QaBO- - Qi^Blfi^^^ " QpBlti^B ' P^) 

B 

-i^^.u;^'^ + \mluJ^u:^' - ^p^, ■ p'^'' + ^mlp^ ■ + £yy. (1) 

The isospin multiplets for baryons B are represented by the Dirac bispinor with vacuum 
baryon mass ms, and isospin operator t^, and a;^,^ and p^u are field strength tensors. The scalar 
self-interaction term [8] 

Uia) = ^g2a' + ^gsa\ (2) 

is introduced to reproduce the correct compressibility of nuclear matter. We perform this calcula- 
tion in the mean field approximation ^ . The Lagrangian density for hyperon-hyperon interaction 
is given by 

B 

+ 1 {d^^*d^a* - mfa*') - ^c/^^.r" + ^m^,/)^^/,^ (3) 
The scalar density and baryon number density are 

s _ 2Jb + 1 [^^b m% 2 



nn = — / o.. /o fe dk , (4) 

^ 27r2 Jo {k^ + mfy/^ ' ^ ^ 



UB = (2Jb + 1)^ , (5) 
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where Fermi momentum is kpg, spin is Jb, and isospin projection is I^b- Effective mass and 
chemical potential of baryon B are = ruB — gaB<7 — 9(7*bO'* and hb = {kpg + m^)^/^ + 
QujB^o + 9<f>B4'o + hBgpBP03, respectively. Charge neutrality in the hadronic phase is imposed 
through the condition 

Q = Y1 QBriB - We - = , (6) 

B 

where is the number density of baryon B, qb is the electric charge and rig and arc charge 
densities of electrons and muons respectively. The total energy density in the hadronic phase is 
given by 

+ E ^ r\k' + mf)'/'edk, (7) 

, _ _ Jo 



and the pressure is 



1 2 2 1 3 -'^ 4 

"2"^^^ - ^920- - -gscr 

l 2 *2,1 22,-'- 2/2,-'- 22 

1 ^ 2 Jb + 1 /■'=^s k"^ dk 
^3 ^ 27r2 /o (fc2 + mf)V2 

^3 , ^ ^ io (A;2 + m2)V2 " 

l=e~,iJ.~ ^ ' ' 



2.2 Antikaon condensed phase 

The constituents of the pure antikaon condensed phase are baryons (neutrons, protons), Icptons 
(electrons, muons) and antikaons, where the baryons are embedded in the condensate. The 
(anti)kaon-(anti)kaon interaction in the pure condensed phase is described using the relativistic 
field theoretical approach, through the exchange of a, u, p, a* and (j) mesons. However, nucleons 
do not couple with the strange mesons, hence g^* n = ^^^r = 0. 

The Lagrangian density for (anti)kaons in the minimal coupling scheme is, 

Ck = DlKD^'K - m*^KK , (9) 

where the covariant derivative = 9^ + iguiK^^n + ig<j>K4'fj. + igpK^K • and the effective mass 
of (anti)kaons is rn^ = nix — QaK^^ — g<j*K(^* ■ The in-medium energies of K~ mesons for s-wave 
(A: = 0) condensation is given by 

pk- = m*K - gtoRi^o - g<pK4>o - \gpKPo3 , (10) 
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where the isospin projection I^^- = —1/2. The scalar and number density of baryon B in the 
antikaon condensed phase are given by 

_ 2J, + 1 f>^^. m*, ^^^^ 



^ 27r2 7o (A;2 + m|2)i/2 

nf = (2Jb + 1)^, (12) 
The scalar density of K~ mesons in the condensate is given by [TO] 

riK- = 2 {^K- + QluKUJo + g,pK4>Q + ij^QpKPo^ KK = 2m*j^KK . (13) 
The total charge density in the antikaon condensed phase is 

= X! IB^B - - ne - = 0. (14) 

B=n,p 

The total energy density in the antikaon condensed phase is 
£ = +^920- +^g3cr +-m^.a + -m^uJo + -m^(t>o + -rUpPos 

B=n,p ^ l=e- ,11- ^ 

+m*j^nK-, (15) 



dk 



where last term denotes the contribution of the K condensate. The pressure in this phase is 

TjK ^22 ^ 3 ^ 4 ^2*2|-'-22|l2i,2|-'^22 

P = -^^^^ -1^920 --gsd -i^rn^,a + -m^uj^ + -m^(l)Q + -rripPQ^ 
1 ^ 2Jb + 1 f^i'B fc^ dk 1 V — f^'^' 



2.3 The Mixed Phase 

The mixed phase of hadronic and condensed matter is governed by the Gibbs phase 
equilibrium rules [TOl E] , 

ph ^ pK^ ^^^^ 

Pb = Pb- (18) 

where /i^ and ji^ are chemical potentials of baryons B in the pure hadronic and K~ condensed 
phase, respectively. The global charge neutrality and baryon number conservation laws are 

(i-x)QVxQ'' = o, (19) 

nB = {l- xh'h + xnf , (20) 

where x is the volume fraction of condensed phase in the mixed phase. The total energy 
density in the mixed phase is given by 

e = (l-x)e" + xe'' • (21) 
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2.4 Bulk Viscosity 

Energy dissipation due to pressure and density variations associated with r-mode oscillation, 
which drive the system out of f3 equilibrium, gives rise to bulk viscosity. The reactions between 
different constituent particles try to bring the system back to an equilibrium configuration, with a 
delay which depends on the characteristic timescale of the interaction. Strong interaction processes 
are insignificant as the strong interaction equilibrium is reached so fast that these processes are 
considered to be in equilibrium compared to typical pulsation timescales. 

As we are concerned about bulk viscosity coefficient in young neutron stars where temperature 
is ~ 10^ - lO^'' K, we want to find out whether non-leptonic processes involving antikaons might 
lead to a high value for the bulk viscosity coefficient. The relevant non-leptonic process involving 
antikaons is 

71^ P + K-. (22) 

As this reaction involves variation of neutron number density (n„) due to density perturbation, 
we consider neutron fraction as a primary parameter. The general expression for the real part of 
bulk viscosity coefficient [Sj [12] is 

1 + [UJT) 

where P is the pressure, r is the net microscopic relaxation time and 700 and 70 are 'infinite' and 
'zero' frequency adiabatic indices respectively. The factor 

nl dP dxn 

can be determined from the EoS. Here Xn = ^ gives the neutron fraction in the equilibrium 
state and = J2b''^b is the total baryon density. In the co-rotating frame, the generic equation 
relating the angular velocity (lo) of {l,m) r-mode to the angular velocity of rotation of the star 

2m ( 



The relaxation time (r) for the process is given by [H] 

1 



6fi 6n^ 



(25) 



Here Sn^ = — is the departure of neutron fraction from its thermodynamic equilibrium 
value in the condensed phase. The reaction rate per unit volume is [39] 

< \Mk\' > kj^J^i 

= — le^rW ' ^ ^ 

where kp^ is the Fermi momentum for neutrons in the condensed phase and < |M/^|^ > is the 
squared matrix element, averaged over initial spins and summed over final spins. The in- medium 
energy of mesons in the condensate is = fi^ = fij^- . If the quantity is calculated 
numerically, as soon as we know the relaxation time, we can calculate the bulk viscosity coefficient. 
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2.5 Matrix element 



Next, we focus on the evaluation of the matrix element for the non-leptonic process (22). In 
general, the matrix element for the decay of a baryon to another baryon and a 0~ meson 
can be written as 

M = u{k2){A + Bj5)u{ki) (27) 

where u{ki) and u{k2) are the spinors of neutrons and protons respectively. Here A is the parity- 
violating amplitude, and B is the parity-conserving amplitude. The squared and spin averaged 
matrix element is given by 

< |7W|2 >= 2pi • k2 + m*^m*p)\Af + {k^ ■ k2 - m*^m*p)\B\'^]. (28) 

For s-wave K~ condensation, k^ =0 =^ \ki\ = \k2\. As fermion momenta lie close to the Fermi 
surfaces, ki-k2 = E1E2 — | A;i 1 1 A;2 1 cos 9 = ^nfJ-p — kF„kFp- This leads to the squared matrix element 

< |7W2| >^ 2[(/i„/ip - kpjF, + m>;)|A|2 + (/x„/ip - kF^kpp - (29) 

Here we apply the weak SU(3) symmetry to the non-leptonic weak decay amplitudes for the 
process (22). The weak decays of the octet hyperons can be described by an effective SU(3) 
interaction with a parity violating (A) and parity conserving (B) amplitudes \14:\ I15|. The weak 
operator is proportional to Gell-Mann matrix Ag to ensure hypercharge violation |Ay| = 1 and 
|AI| = 1/2. Similarly, the amplitudes for (22) are extracted from experimentally known de- 
cay parameters of the weak decay of hyperons [15]. The amplitudes are A = —1.62 x 10"'' and 
B = —7.1 X 10~^. It is to be noted that all quantities in (29) are to be calculated in the condensed 
phase. 



2.6 Critical Angular Velocity 

The bulk viscosity damping timescale (tb) due to the non-leptonic process involving antikaons is 
given by [HI [3] 

I 1 dE 

= -TEH • <™> 

where E is the energy of the perturbation as measured in the co-rotating frame of the fluid and 
is expressed as 

E = -a^n^R-^ [ e{r)r^dr . (31) 
2 Jo 

Here, a is the dimensionless amplitude of the r-mode, R is the radius of the star and e(r) is the 
energy density profile. The derivative of the co-rotating frame energy with respect to time is 

_ = -47r / C(^) < |V • 6v\^ > r^dr, (32) 
dt Jo 



where the angle average of the square of the hydrodynamic expansion [T7] is 

690 \RJ \ ^ ' \RJ I \ TrGe ' 



< |V • (Jul^ >= ^— ^ -- 1 + 0.86' ' 
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and e is the mean energy density of a non-rotating star. The total r-mode time scale (r,.) is defined 
as 

1111 

- = + — + — . (33) 

where the time scales for gravitational radiation {tgr) and modified Urea process {tjj) involving 
only nucleons have also been included. The gravitational radiation timescale is given by [18j 



1 1310727r^6 /"^ / N 6j .o.^ 
-ir / e{r)r^dr . (34) 



TGR 164025 Jo 

The time scale due to modified Urea process [tu) involving only nucleons is calculated from (31) 
using the following expression for bulk viscosity coefficient for modified Urea process [181 E] 

= 6 X lO-^^e^T^w^ _ (35) 

For a star of given mass, solving = 0, we can obtain the critical angular velocity at each 
temperature above which the r-mode becomes unstable. 

3 Parameters of the Theory 

3.1 Nucleon-Meson coupling constants 

Nucleon-meson coupling constants are determined from saturation properties of nuclear matter 
[20] . We used the following values: binding energy = —16.3 MeV, baryon density no = 0.153 fm~'^, 
asymmetry energy coefficient Oasy = 32.5 MeV, incompressibility K = 240 MeV, and effective 
nucleon mass m*j^/mN = 0.78. We also studied the parameter set [21] with incompressibility 
K = 300 MeV, and effective nucleon mass m'^/mN = 0.70. 

3.2 Kaon-Meson coupling constants 

According to the quark model and isospin counting rule, the vector coupling constants are given 
by ^ 

9ioK = -^guN and gpx = QpN ■ (36) 

The scalar coupling constant is obtained from the real part of optical potential depth at 
normal nuclear matter density 

i'^o) = -9aKO- - QujK^O ■ (37) 

It is known from /('"-atomic data that antikaons experience an attractive potential in nuclear 
matter while kaons feel a repulsive interaction [22\ [23l [M] [25] [26l [27] . The strength of antikaon 
optical potential depth ranges from shallow attractive (—40 MeV) to strongly attractive (—180 
MeV). Here we perform the calculation for antikaon optical potential depth at normal nuclear 
matter density C/^(no) = —120 MeV. The strange meson fields couple with (anti)kaons. The 
a*-K coupling constant is ga*K = 2.65 as determined from the decay of /o(925) meson and the 
vector (p meson coupling with (anti)kaons V^g(j,K = 6.04 follows from the SU(3) relation [28] . 
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3.3 Hyperon-Meson coupling constants 

Hyperon-vector meson coupling constants are determined from SU(6) symmetry of the quark 
model [291 EQl [31]. The scalar a meson coupling to hyperons is calculated from the potential 
depth of a hyperon (Y) 

up (no) = -gaYcr + 5<^yu;o , (38) 

in normal nuclear matter. The potential depth of A hyperons in normal nuclear matter C/^(no) = 
—30 MeV is obtained from the analysis of energy levels of A hypernuclei [30^ I32j . Recent H- 
hypernuclei data from various experiments [33\ [M] give a relativistic potential of Uff (no) = — 18 
MeV. However, the analysis of S~ atomic data implies a strong isoscalar repulsion for S~ hyperons 
in nuclear matter [35]. Also, recent S hypernuclei data indicate a repulsive E-nucleus potential 
depth [36]. Therefore, a repulsive potential depth of 30 MeV for S hyperons [35] is adapted. 

(y') 

The hyperon-(T* coupling constants are estimated by fitting them to a potential depth, Uy (no), 
for a hyperon (Y) in a hyperon (Y ) matter at normal nuclear matter density obtained from double 
A hypernuclei data [371 [29]. This is given by 

ui\no) = uPino) = 2ui^Hn^) = 2[/f ^(no) = -40 . (39) 

4 Results and Discussions 

The Equation of State (pressure versus energy density) for neutron star matter with K~ con- 
densate is plotted in Fig. 1 (solid line) for K=240 MeV. The equations of state for neutron star 
matter containing nucleons only (short dashed line) and with hyperons (long dashed line) are 
superimposed on the same figure. The two kinks on the EoS involving the condensate at 3.26 no 
and 4.62 no mark the beginning and end of the mixed phase. For K~ condensed matter with 
K=300 MeV, similar kinks are observed at 2.23 no and 3.59 no defining the mixed phase. The 
threshold density for the appearance of the A hyperon is 2.6 no- The EoS becomes softer in 
presence of exotic matter {K^ condensate or hyperons) compared with that of nucleon matter. 
However, the EoS for antikaon condensed matter is stiffer than that for hyperon matter beyond 
the mixed phase. 

For the calculation of damping time scale due to bulk viscosity using (30), we need the energy 
density profile and bulk viscosity profile of the neutron star. We choose a neutron star of gravita- 
tional mass 1.63M0 corresponding to a central baryon density 3.94 no and rotating at an angular 
velocity $1 = 1180s~^. The neutron star is so chosen to ensure that it contains K~ condensate 
in its core because the central baryon density is well above the threshold of K~ condensation 
(3.26 no). The bulk viscosity profiles for K = 240 MeV and K = 300 MeV are plotted in Fig 2 
as a function of equatorial distance. Here we note that the bulk viscosity profile drops to zero 
value beyond a certain equatorial distance, when the baryon density decreases below the threshold 
density of condensation and the non-leptonic process in (22) ceases to occur in the star. 

We compare the bulk viscosity coefficient due to antikaon condensate with bulk viscosity due 
to the non-leptonic weak interaction n+p ^ p + K involving hyperons. The coefficient of bulk vis- 
cosity due to non-leptonic process involving hyperons is shown in Fig 3 as a function of normalised 
baryon density for different temperatures [38j. Comparing the two figures 2 and 3, we can infer 
that the bulk viscosity coefficient in antikaon condensed matter is suppressed by several orders of 
magnitude in comparison with nonsuperfluid hyperon bulk viscosity. It must also be noted that 
the antikaon bulk viscosity is independent of temperature, while hyperon bulk viscosity increases 
with decrease in temperature. Hence, hyperon bulk viscosity can act as an effective damping 
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mechanism at low temperatures. 

We can obtain critical angular velocities as a function of temperature by solving = in (33) 
for a rotating neutron star of mass 1.63 M solar- We calculate the bulk viscosity coefficient due 
to modified Urea process due to nucleons using (35). The modified Urea bulk viscosity is plotted 
as a function of normalised baryon density for a range of temperatures in Fig 5. It is evident 
from the figure that bulk viscosity due to modified Urea process involving nucleons increases with 
increase in temperature. Hence r-mode instability is effectively suppressed by modified Urea bulk 
viscosity at high temperatures. The ratio of critical angular velocities to the rotational velocity 
of the neutron star are plotted as a function of temperature in Fig 4 for non-leptonic processes 
involving antikaons as well as hyperons. In this figure, we observe that there exists a window of 
instability for matter with hyperons. We can interpret that hyperon bulk viscosity damps the 
r-mode instability at low temperatures, whereas it is effectively suppressed at high temperatures 
by the bulk viscosity due to the modified Urea process involving nucleons. For K~ condensed 
matter, we can infer that r-mode instability is damped by nucleonic modified Urea process at 
high temperatures. From the figure, it is evident that antikaon bulk viscosity is not an effective 
mechanism to damp the instability at low temperatures. The argument is further justified by the 
fact that the critical velocity curves obtained for both parameter sets K= 240 MeV and K = 300 
MeV are dictated purely by the bulk viscosity due to nucleonic modified Urea process. The bulk 
viscosity due to antikaon condensation is not sufficient to suppress the r-mode within the tem- 
perature range considered here. Hence the instability window is wider for K~ condensed matter, 
and neutron stars with K~ condensate in its interior could be possible sources of gravitational 
waves. 

5 Summary 

The role of K~ condensation on bulk viscosity and r-mode instability has been investigated in 
this paper. We have estimated the bulk viscosity coefficient and the corresponding damping time 
scale due to the non-leptonic process n ^ p + K~ and compared them with those associated with 
the non-leptonic process n + p ^ p -\- A. involving hyperons. We have considered a first order 
phase transition from the nuclear to the antikaon condensed phase, and the equation of state has 
been constructed using relativistic mean field theoretical models. We find that the bulk viscosity 
coefficient in K~ condensed is suppressed by several orders of magnitude in comparison with the 
nonsuperfluid hyperon bulk viscosity [3l [T61. 138j . We also infer that the antikaon bulk viscosity 
is unable to damp the r-mode instability unlike hyperon bulk viscosity. If the instability is not 
suppressed by the viscosity, the star is forced to lose angular momentum as it cools down via 
gravitational radiation. So, r-modes in neutron stars with antikaon condensate in their core are 
possible candidates of gravitational waves. 
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Fig. 1. The equation of state ( pressure P vs energy density e ) for matter containing nuclcons only (short 
dashed hue), with hyperons (long dashed line) and with K~ condensate (solid line) for antikaon optical 
potential depth at normal nuclear matter density Uj^ = —120 MeV. 
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Fig. 2. Bulk viscosity profile is plotted with equatorial distance for a rotating neutron star of mass 1.63 
M0 for two parameter sets. 
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Fig. 3. Hyperon bulk viscosity coefBcient is exhibited as a function of normalised baryon density for 
different temperatures. 



14 




Fig. 4. Density dependence of modified Urea bulk viscosity is shown for a range of temperatures. 
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Fig. 5. Critical angular velocities for 1.63 Mq neutron star are plotted as a function of temperature. 
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